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Abstract 

, Bianchi type V bulk viscous fluid cosmological models are investi- 

U" 1 gated with dynamic cosmological term A(t). Using a generation tech- 

, 1 | ' nique (Camci et ai, 2001), it is shown that the Einstein's field equations 

(3JT), are solvable for any arbitrary cosmic scale function. Solutions for partic- 

JL" ■ ular forms of cosmic scale functions are also obtained. The cosmological 

. ' constant is found to be decreasing function of time, which is supported by 

results from recent type la supernovae observations. Some physical and 

5_l ' geometrical aspects of the models are also discussed. 
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1 Introduction 

The study of Bianchi type V cosmological models create more interest as these 
models contain isotropic special cases and permit arbitrary small anisotropy lev- 
els at some instant of cosmic time. This property makes them suitable as model 
of our universe. The homogeneous and isotropic Friedman-Robertson- Walker 
(FRW) cosmological models, which are used to describe standard cosmological 
models, are particular case of Bianchi type I, V and IX universes, according 
to whether the constant curvature of the physical three-space, t — constant, is 
zero, negative or positive. These models will be interesting to construct cos- 
mological models of the types which are of class one. Present cosmology is 
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based on the FRW model which is completely homogeneous and isotropic. This 
is in agreement with observational data about the large scale structure of the 
universe. However, although homogeneous but anisotropic models are more re- 
stricted than the inhomogeneous models, they explain a number of observed 
phenomena quite satisfactorily. This stimulates the research for obtaining exact 
anisotropic solution for Einstein's field equations (EFEs) as a cosmologically 
accepted physical models for the universe (at least in the early stages). Roy 
and Prasad [T] have investigated Bianchi type V universes which are locally 
rotationally symmetric and are of embedding class one filled with perfect fluid 
with heat conduction and radiation. Bianchi type V cosmological models have 
been studied by other researchers (Farnsworth [2J , Maartens and Nel [5J , Wain- 
wright et al. 4 , Collins [5], Meena and Bali [6J, Pradhan et al. [3 Ej) in 
different context. 

Models with a dynamic cosmological term A(t) are becoming popular as they 
solve the cosmological-constant problem in a natural way. There is a significant 
observational evidence for the detection of Einstein's cosmological constant, A or 
a component of material content of the universe, that varies slowly with time and 
space and so acts like A. Recent cosmological observations by High-z Supernova 
Team and Supernova Cosmological Project (Garnavich et al. [OJ, Perlmutter et 
al. [TO], Riess et al. [TT], Schmidt et al. [12]) strongly favour a significant 
and positive A with the magnitude A(Gh/c 3 ) ss 10~ 123 . These observations on 
magnitudes and red-shift of type la supernova suggest that our universe may 
be accelerating one with a large function of the cosmological density in the 
form of cosmological A-term. Earlier researches on this topic, are contained 
in Lodovico [15] . Weinberg [14], Dolgov [15] — [17], Bertolami [18], Ratra and 
Peebles [19], Carrol, Press and Turner [20] . Some of the recent discussions on the 
cosmological-constant "problem" and consequence on cosmology with a time- 
varying cosmological-constant have been discussed by Tsagas and Maartens [21] , 
Sahni and Starobinsky [22], Peeble [23], Padmanabhan [24], Vishwakarma [25] , 
and Pradhan et al. 26J. This motivates us to study the cosmological models, 
where A varies with time. 

The distribution of matter can be satisfactorily described by a perfect fluid 
due to the large scale distribution of galaxies in our universe. However, observed 
physical phenomena such as the large entropy per baryon and the remarkable 
degree of isotropy of the cosmic microwave background radiation, suggest anal- 
ysis of dissipative effects in cosmology. Furthermore, there are several processes 
which are expected to give rise to viscous effects. These are the decoupling of 
neutrinos during the radiation era and the decoupling of radiation and matter 
during the recombination era. Bulk viscosity is associated with the GUT phase 
transition and string creation. Misner [27] has studied the effect of viscosity on 
the evolution of cosmological models. The role of viscosity in cosmology has 
been investigated by Weinberg [25]. Nightingale [25], Heller and Klimek [50] 
have obtained a viscous universes without initial singularity. The model stud- 
ied by Murphy |5T| possessed an interesting feature in which big bang type of 
singularity of infinite space-time curvature does not occur to be a finite past. 
However, the relationship assumed by Murphy between the viscosity coefficient 
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and the matter density is not acceptable at large density. Thus, we should con- 
sider the presence of material distribution other than a perfect fluid to obtain a 
realistic cosmological models (see Gr0n [33] for a review on cosmological models 
with bulk viscosity) . The effect of bulk viscosity on the cosmological evolution 
has been investigated by a number of authors in the framework of general theory 
of relativity. This motivates to study cosmological bulk viscous fluid model. 

In recent years, several authors (Hajj-Boutros [33], Hajj-Boutros and Sfeila 
[33], Ram [35], Mazumder [3J5] and Pradhan and Kumar [37]) have investi- 
gated the solutions of EFEs for homogeneous but anisotropic models by using 
some different generation techniques. Bianchi spaces I — IX are useful tools in 
constructing models of spatially homogeneous cosmologies (Ellis and MacCal- 
lum [38], Ryan and Shepley [39]). From these models, homogeneous Bianchi 
type V universes are the natural generalization of the open FRW model which 
eventually isotropize. Recently Camci et al. 40J derived a new technique for 
generating exact solutions of EFEs with perfect fluid for Bianchi type V space- 
time. Very recently Pradhan et al. [H] have obtained Bianchi type V perfect 
fluid cosmological models with time dependent A-term. 

In this paper, in what follows, we will discuss Bianchi type V cosmological 
models obtained by augmenting the energy-momentum tensor of a bulk viscous 
fluid by a term that represents a cosmological constant varying with time, and 
later generalize the solutions of Refs. [35] 001 SI]- This paper is organized 
as follows: The field equations and the generation technique are presented in 
Section 2. We relate three of the metric variables by solving the off-diagonal 
component of EFEs, and find a second integral which is used to relate the 
remaining two metric variables. In Section 3, for the particular form of each 
metric variables, some solutions are presented separately and solutions of Camci 
et al. [20], Ram [35] and Pradhan et al. [3J are shown to be particular cases 
of our solutions. Kinematical and dynamical properties of all solutions are also 
studied in this section. In Section 4, we give the concluding remarks. 



2 Field equations and generation technique 

In this section, we review the solutions obtained by Pradhan et al [H] . The 
usual energy-momentum tensor is modified by addition of a term 

T^ = -k(t) 9l3 , (1) 

where A(i) is the cosmological term and gij is the metric tensor. Thus the new 
stress energy-momentum tensor in presence of bulk stress is given by 

Tij =(p + p)uiUj - pgij - A(t)gij, (2) 

where 

P = P + H (3) 
Here, p, p, p, £ and u are, respectively, the energy density, isotropic pressure, 
effective pressure, bulk viscous coefficient and the fluid four-velocity vector of 
distribution such that u l Ui = 1. In general, £ is a function of time. 
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We consider the space-time metric of the spatially homogeneous Bianchi 
type V of the form 



ds 2 = dt 2 - A 2 (t)dx 2 



JlOLX 



B 2 {t)dy 2 + C 2 (t)dz 2 



(4) 



where a is a constant. For the energy momentum tensor ^ and Bianchi type 
V space-time @, Einstein's field equations 



yield the following five independent equations 
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Here and in what follows the suffix 4 by the symbols A, B, C and p denote 

form 

Pi + (p + P )e = o. (11) 



differentiation with respect to t. The Bianchi identity (T?j — 0) takes the form 



It is worth noting here that our approach suffers from a lack of Lagrangian 
approach. There is no known way to present a consistent Lagrangian model 
satisfying the necessary conditions discussed in this paper. 

The physical quantities expansion scalar 6 and shear scalar a 2 have the 
following expressions: 

> tt i = ^ + ^ + S (12) 



1 a 1 

-ana 3 = - 

2 3 3 
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Integrating Eq. 
obtain 



AB AC BC 
TO]) and absorbing the integration constant into B or C, we 



A = BC 



without any loss of generality. Thus, elimination of p from Eqs. 
the condition of isotropy of pressures 
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which on integration yields 

£>4 C4 



B c {Bey/ 2 ' 



(16) 



where k is a constant of integration. Hence for the metric function B or C from 
the above first order differential Eq. (flB|) , some scale transformations permit us 
to obtain new metric function B or C. 

Firstly, under the scale transformation dt = B x / 2 dT, Eq. (p~6|) takes the form 

CB T - BC T = kC~ 1/2 , (17) 

where subscript represents derivative with respect to r. Considering Eq. (JTTJ) as 
a linear differential equation for B, where C is an arbitrary function, we obtain 

f dr 

(f)B = k 1 C + kCj -g^, (18) 

where fej is an integrating constant. Similarly, using the transformations dt = 
B 3 / 2 df, dt = C l l 2 dT, and dt = C 3 ^ 2 df in Eq. (Ig), after some algebra we 
obtain respectively 

(it)fl(r; fc 2 , fc) - fc 2 C exp f * | , (19) 

/dT 
&J2> ( 20 ) 



and 



(iv)C(f; h, k) = hB cxplkj -J^ j , (21) 



where &2, &3 and &4 are constants of integration. Thus choosing any given 
function B or C in cases (i), (ii), (iii) and (iv), one can obtain B or C and hence 
A from IE 



3 Generation of new solutions 

We consider the following four cases: 

3.1 Case (I): Let C = r n , (n is a real number satisfying 

In this case, Eq. (fl8j) gives 

B = hr n + J^ r i-3»/2 (22) 
2-5?i v 7 
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and then from (fT4"]) . we obtain 

A 2 = k lT 2n 



2k 



-l-n/2 



(23) 



2 - 5n 

Hence the metric ((H) reduces to the new form 

ds 2 = (fcir™ + 2£r ei ) [dr 2 - T n dx 2 } - e 2ax Uk lT n + 2ir^) 2 dy 2 + T 2n dz : 
where 

l = 2 J ^ andh = l - i i • 

For this derived model (|24p . the effective pressure, energy density and cosmo- 
logical constant are given by 



(24) 



8tt(p + A) = (feir" + 2£r ei ) 



-2fefn(n - l)r 2 "~ 2 - Mn(lG - 13n)T-( fl+2 ") 



^ 2 (4 + 4n- lln 2 ) ^_ 3 „ 



, -3 



8tt(p- A) = (fc 1 r n +2^r fl 

£ 2 (4 + 4ri- lln 2 ) _ 3n 



+ a 2 r~" (jfciT" + 2£t £i ) , (25) 
3fc 2 n 2 r 2 "- 2 + 3fciM2 - n)T-^ +2n ^ 



Za 2 T- n (fcir™ + 2£t £ 



(26) 



We assume for the specification of £, that the fluid obeys an equation of state 
of the form 

P = IP, (27) 

where 7(0 < 7 < 1) is a constant. 

Thus, given £(i), we can solve the system for the physical quantities. In 
most of investigations involving bulk viscosity is assumed to be a simple power 
function of the energy density (Pavon et aZ. [42] ; Maartens [43]; Zimdahl [44] ; 
Santos et al. |45j ) 

m=tioP w , (28) 

where £0 and w are real constants. For small density, w may even be equal to 
unity as used in Murphy's work 46J for simplicity. If w — 1, Eq. ([28)1 may 
correspond to a radiative fluid (Weinberge [H]). Near the big bang, < w < \ 
is a more appropriate assumption (Belinskii and Khalatnikov [47]) to obtain 
realistic models. 

For simplicity and realistic models of physical importance, we consider the 
following two cases (w = 0, 1): 
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3.1.1 Model I: Solution for w = 

When w = 0, Eq. I|28p reduces to £ = £o = constant. Hence in this case Eq. 
Z3J), with the use of J27]) and ([25]). leads to 



8?r(l + 7)p = Z?r 3 L(n + 2)fc 2 T 2(n ~ 1) + 2(5n - 2)M«T^ 1+2n) 



2a 2 r- n + 24tt£o ( /fcinr™^ 1 + \i{2 - njT'^ ) D x 



where 



(29) 



Di =k 1 T n + 2£r ei . 
Eliminating p(t) from Eqs. (f26|) and (|29| . we obtain 

8vr(l+7)A = L>r 3 [(l-2n-3n7)nfc2 T 2 (»- 1 ) -{(10-13n)+3(2-n)7}fcifo"-- (£l+2n) 



-_(4 + 4n - lln 2 )(l + 7^ 2 r- 3 " 



(1 + 37)a 2 T-"L>r 1 + 247r$o 



n/cir"" 1 + -1(2 - n)T 



(30) 



3.1.2 Model II: Solution for w = 1 

When w — 1, Eq. (j2"5j) reduces to £ = £oP- Hence in this case Eq. (|2"5|) . with 
the use of (HI, ([27]) and JUL leads to 



8?rp 
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(31) 



Eliminating p(i) between ([2"o| and ([3"Tj) , we obtain 
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(32) 

From Eqs. (f2"9"| and ([3"0")l . we observe that at the time of early universe the 
the energy density p(t) and cosmological constant (A(i)) decrease with time in- 
creases (see Figure 1). We also observe that the value of A is small and positive 
at late times which is supported by recent type la supernovae observations [2j 
— [T2]. In Model II, from Eqs. ([3T[) and ([3"2"|) . we find that for large range of 
parameters the energy density decreases with time very sharply and becomes 
negative (even if it is positive initially) and then remains negative throughout 
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the evolution. The cosmological constant shows singular behaviour closer to ori- 
gin. Then it shows eretic behaviour and later stage it remain negative constant 
value. It seems that Model II may not be physical model of the universe. 

The metric (f2"4"]l is a four-parameter family of solutions to EFEs with a bulk 
viscous fluid. Using the scale transformation dt = B^dr in Eqs. (fl~2]) and (TT3"]) 
for this case, the scalar expansion 6 and the shear a have the expressions: 



klTLT 
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Eqs. and (O lead to 



a k 
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(33) 
(34) 

(35) 



Now, we consider four subcases for the parameters A, n, k, k\ whether zero or 
not. 



In subcase A = and £o = 0, metric ([24]) with expressions p, p, 6 and a 
for this model are same as that of solution (fl~8|) of Camci et al. [40] . If we set 
£o = 0, the metric ([24]) gives the solution obtained by Pradhan et al. [41] . 

In subcase A = 0, n — 0, after a suitable inverse time transformation, we 
find that 

ds 2 = dt 2 - K x {t + t^f^dx 2 - e 2ax \Ki(t + hf^dy 2 + dz 2 ] , (36) 

where to is a constant of integration and K\ = (3fc/2) 2 / 3 . The expressions p, p, 
9 and a for this model are not given here, since it is observed that the physical 
properties of this one are same as that of the solution ([24]) of Ram [35 . 



In subcase A = 0, k = 0, after inverse time transformation and rescaling, 
the metric (f2"4"]l reduces to 

ds 2 = dt 2 - K 2 (t + ti)?& [dx 2 + e 2ax (dy 2 + dz 2 )] , (37) 

where t\ is a constant of integration and K 2 — (^j ) " +2- ^ or som ti° n 5 
when n = 1 and a = 0, we obtain Einstein and de Sitter [48] dust filled universe. 
For K 2 = 1, t\ = and n = ' wnere w is a parameter in Ram's paper 

[35], the solution ([37]) reduces to the metric (|14l) of Ram [35]. In later case, if 
also a — 0, then we get the Minkowski space-time. 
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Now, in subcase A = 0, fci = 0, after some algebra the metric (|2"4"]) takes the 
form 

ds 2 = dt 2 - 2£K 3 {t + t 2 ) 2 / 3 [dx 2 + e 2ax (at mi dy 2 + a^r^dz 2 )] , (38) 

where t 2 is a constant, at mi = 2£K^ 7r (t + t 2 )^=$ and K 3 = 3{2 ^ 23 . 
For t% = 0, k = I and n = from (|3"5|) . we obtain that the solution (J2H) of Ram 

1351 • 

Some Physical aspects of model : 

The model (|24|) has barrel singularity at r = tq given by 



(2-5n) 



which corresponds to t = 0. For n ^ 2/5 from it is observed that at 

the singularity state r = To, p, p, A, and a are infinitely large. At t — > oo, 
which corresponds to r — > oo for n < 2/5 and > 0, or r — » for n > 2/5 
and k < 0, p, p, A, and cr vanish. Therefore, for n ^ 2/5, the solution (f24|l 
represents an anisotropic universe exploding from r = to, i.e. < = 0, which 
expands for < t < oo. We also find that the ratio er/(9 tends to a finite limit 
as t — > oo, which means that the shear scalar does not tend to zero faster than 
the expansion. Hence the model does not approach isotropy for large values of 
t. 

In subcase A = 0, k = 0, the ratio (|35| tends to zero, then the model 
approaches isotropy i.e. shear scalar a goes to zero. For the model (|3"T1) . p and 
p tends to zero as t — > oo; the model would give an essentially empty universe 
at large time. The dominant energy condition given by Hawking and Ellis [49] 
requires that 

p + p>0, p + 3p>0 (39) 

Thus, we find for the model ([3T]) that n(2 ~ n) > 0. Hence for the values 
< n < 2, the universe |37|) satisfies the strong energy condition i.e. p + 3p > 0. 
Also this model is sheer-free and expanding. 

In subcase A = 0, k\ = 0, for n ^ 2/5, 2, it is observed from relations (I29|) 
- (|34p that p, p, 9 and u are infinitely large at the singularity state t = —t 2 . 
When t — > oo, these quantities vanish. We also find that the ratio a/9 is a con- 
stant. This shows that the cosmological model ([38)) does not approach isotropy 
for large value of t. In this model the dominant energy conditions ()39|) are then 
verified for 6 — 5n — 25n 2 > 0. Since n ^ 2/5, the model ([38| satisfies the strong 
energy condition for —3/5 < n < 2/5. 

In each of subcases, all the obtained solutions (f36|) . ([37)) and (|38|) satisfy the 
Bianchi identity given in Eq. (|10p . 



to 



fci(5n- 2) 
2k 
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3.2 Case (II): Let C = f n , (n is a real number satisfying 

n^2/3). 

In this case Eq. (fT9|) gives 

B = k 2 T n exp (Mf £l ) (40) 

and from (fT4"]) , we obtain 

A 2 = fc 2 f 2 " exp (Mf £l ) (41) 
where M = j-. Hence the metric (0} reduces to the form 



ds 2 = f^-W* 



f 2{l-l,)/3 (? MT^ d ~2_ e MT*r dx 2 



e 2ax (e 2Mftl dy 



2 ' dz 2 



(42) 



where the constant k 2 is taken, without any loss of generality, equal to 1. This 
metric is a three-parameter family of solutions to EFEs with a perfect fluid. 

For the above model, the distribution of matter and nonzero kinematical 
parameters are obtained as 

Mp + A) = 2nf 2 ^-^ + 3nkf 3e ^ 

+ £fHh-l) + Q 2 f 4(fe-l)/3 e -3Mf^ _ (43) 



3^1-4 



8vr(p - A) = 3n 2 f 2 ^^ + 3nkf 
+ ^f^-V - ^fMtx-iySg-SMf*! (44) 

For simplicity and realistic models of physical importance, we consider the 
following two cases (w = 0, 1): 

3.2.1 Model I: Solution for w = 

When w = 0, Eq. I|28p reduces to £ = £o = constant. Hence in this case Eq. 
1(351). with thc usc of 62)' <G3 and (ESI), leads to 

8tt(1+7)p= n(2 + 3n)f 2(£l ~ 2) +6nfcf (3£l " 4) + fc 2 f 4(£l " 1) -2L)3 + 24<o£ ) 2, (45) 
where 

D 2 = nf^-V + ifef2(A-l) ) 

Z?3 = a 2 f^ £l - 1 )e- 3Mf ' 1 . 
Eliminating p(i) from Eqs. ([44]) and (|45| . we obtain 

8tt(1 + 7 )A = n(2 - 3n 7 )f 2( ^ 2) + 3nfc(l - 7 )f (3£l ~ 4) 

+ ifc 2 (l - 7 )f 4 (^-D + (1 - 3 7 )£ 3 + 24^oA>. (46) 
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3.2.2 Model II: 



Solution for w = 1 



When w = 1, Eq. (|28[) reduces to £ = £oP ■ Hence in this case Eq. (j4"3"]) . with 
the use of (gl|), {27]) and (gSJ), leads to 

n(2 + 3n)f 2 ^" 2 ) + 6nfcf ( 3£ ^ 4 ) + k 2 f 4 ^~V -2D 3 
87TP = (1 + 7)-%^ ' (4?) 

Eliminating p(t) from Eqs. f4"4"]) and l(4"7]) , we obtain 

n(2 + 3n)f 2 ^- 2 ) +6nfcf( 3 ^- 4 ) +fc 2 f 4 ^- 1 ) -2Z?3 

o7rA = — — 

(l+7)-3^£>2 

- in^2)f 2 ^-^ - Snkf^-V - l -k 2 r^-^ + 3D 3 . (48) 

In Mode I, from Eqs. (f4"5]) and (HBJ), we observe that the energy density p(t) 
and cosmological constant (A(t)) decrease with time increases (see Figure 2). 
Here we find energy density always positive. We also observe that the value 
of A is small and positive at late times which is supported by recent type la 
supernovae observations [9] — [12]. From Eqs. (|47|) and 1(35]). we observe that 
Model II has the similar behaviour as Model I, so it is not reproduced here. 

The scalar of expansion 9 and the shear a are obtained as 

k. 



= 3 

fc 
2 

From Eqs. P9"]) and lp)0]) . we have 



n f*i-* + If 2(^-1) 



(49) 



CT =|f 2 ^- 1 ) e - 3Mfl , (50) 



6(nf-^ + |; 



(51) 



In subcase A = and £o = 0, metric ((42]) with expressions p, p, 9 and cr for 
this model are same as that of solution ((27]) of Camci et al. |40j . 

In sub-case A = 0, £o = 0, l\ = 1 (i.e. n = 0), we find a similar solution to 
((31)]) , and hence this subclass is omitted. For k = 0, the ratio (jST]) is zero and 
hence there is no anisotropy. 



After a suitable coordinate transformation, the metric Q42J can be written 

as 

ds 2 = dt 2 - K 4 (t + t 3 ) 2ei [dx 2 + e 2ax (dy 2 + dz 2 )] , (52) 
where t 3 is a constant and K4 



Mi = 2^ 7^ |, where Mi is a 



(2-3Mi) 

new parameter. When M% =0 and l\ = 6, from ((52"]) , we get the solution |T 
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of Ram [35] . 



Some physical aspects of model: 

The models have singularity at f — > — oo for £\ > or f — > for £\ < 0, 
which corresponds to t — > 0. It is a point type singularity for £\ > whereas it is 
a cigar or a barrel singularity according as £\ < 0. At t — > oo, which correspond 
to f — > oo for 4 > or f — > for 4 < 0, from Eqs. ([45]) - ([50]), we obtain that 
for £i > 0, p, p — > 0, and er, 6> — > (A; > 0), - oo (fc < 0); for £ 1 < 0, similar the 
above ones. Then, clearly for a realistic universe, it must be fulfill as f — > — oo, 
n and k are positive and i\ is an odd positive number; as f — > 0, fc is positive, 
and t\ an even negative number. Also, since lim^—,^^ ^ 0, therefore these 
models do not approach isotropy for large values of f . 



In subcase k — for the metric ([52]) . the effective pressure, density and 
cosmological constant are given by 

AN £i(2-3£i) a 2 
3^i 3q 2 

8 " (p - A)= (tT^- ^(t + t 3 )" 1 - (54) 

When A = and ^ = 0, the pressure and energy density are same as that 
of given in Eq. ([4~4"]) of paper Camci et al. [ID] . In this case, the weak and 
strong energy conditions ([39]) for this solution are identically satisfied when 
^1(1 — l\) > i.e. < £1 < 1. This model is shear-free and expanding with 

- 3£l 

3.3 Case (III) : Let B = T n (n is a real number). 

In this case Eq. ([20]) gives 

C = k 3 T n - 2£T ei (55) 

and then from (Q3J), we obtain 

A 2 = k 3 T 2n - 2£T ll+n (56) 
Hence the metric ((U) takes the new form 

ds 2 = (k 3 T n - 2£T £l ) [dt 2 ~ T n dx 2 }- 

T 2n dy 2 + (k 3 T n - 2£T l1 ) 2 dz 2 



(57) 



For four-parameters family of solution (|57[) . the physical and kinematical quan- 
tities are given by 



8ir(p + A) 



£ 2 

— (lln 2 -An- 4)T- 3n + £k 3 n(13n - 10)T fl+n - 
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2k 2 3 n{n - \)T 2 



8ir(p - A) = 



3kln 2 T 2n ~ 2 



(k 3 T n - 2£T £l ) 
l 2 {\\n 2 -An- 14) 



2 T- n (k 3 T n - 2£T^) 



(58) 



T~ 4n - 3£k 3 n(2 - n)T 



ti+n 



(k 3 T n - 2£T gl ) 3 - 3a 2 T~ n (k 3 T n - 2lT ll \ 



(59) 



When w = 0, Eq. (|28[) reduces to £ = £o = constant. Hence in this case Eq. 
(|58j) , with the use of ([59]), J23 and dHJ) , leads to 



3.3.1 Model I: Solution for w = 



8tt(1+7)p = \n(n + 2)k 3 2 T 2{n - l) + 16n(n - l)k 3 £T^ 1+n) - {lln 2 -An- 9)£ 2 T~ 

(60) 



where 



2a 2 T~ n D 4 r 1 + 247r£ A>£>4~ 5 . 
L> 4 = k 3 T n - 2£T £l , 



D 5 = fonT 71 ' 1 + ^{n- 2)£T~^ . 
Eliminating p(t) from Eqs. ([59")) and (fBT)]) . we obtain 

8tt(1 + 7)A = {-Ki£ 2 {\ - 7)(lln - 4) - ^(77 - 2)}T" 



3n 



+n£fc 3 {13n - 10 + 3 7 (2 - n)}T ll+n - nk 3 2 (2n -2 + 3n 7 )T 2 ("- 1) 
+ (1 + 3j)a 2 T- n D 4 - n + 247r£oA>AT=- 

3.3.2 Model II: Solution for w = 1 



£4 



-3 



(61) 



When m; = 1, Eq. ([28]) reduces to £ = £oP- Hence in this case Eq. ([58]) . with 
the use of (J5SJ), ([271) and (gSJ), leads to 



8?rp 



£ 2 (lln 2 - 4n - 9)T" 3n + 16n(n - l)lk 3 T ll+n 

(62) 



[(1 + 7)-3^ 5j D4" 5 ] 

+ n(n + 2)k 3 2 T 2{n - 1) ] L> 4 -3 - c^T^AT 1 . 
Eliminating p(t) from Eqs. ([59"]) and (|6"2)) , we obtain 
1 

[(l + 7 )-3£ IW 

+n(n + 2)/c 3 2 T 2 ("- 1 )}L> 4 ~ 3 - a 2 T~ n D 3 ^ 1 



8ttA = 



{-^ 2 (lln 2 -An- 9)T~ 3n + 16n(n - l)£k 3 T ll+n 
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{-i(lln 2 -An- U)£ 2 T~ 3n - 3n(2 - n)£k 3 T ll+n + 3n 2 k 3 2 T^ n -^}Dr 3 + 3a 2 T"L> 4 _1 

(63) 

In Mode I, from Eqs. (fBT)]) and (|6T|) . we observe that the energy density p(t) 
and cosmological constant (A(t)) are decreasing functions of time (see Figure 
2). The energy density is always positive. We also observe that the value of A is 
small and positive at late times which is supported by recent type la supernovae 
observations [9] — [12]. From Eqs. (|62|) and (|63| . we observe that Model II has 
the similar behaviour as Model I, so it is not reproduced here. 



The scale of expansion and the shear are obtained as 



= 3 



£{n 2) T' 3 "/ 2 + k 3 nT n 
2 J 



l.rp-3n/2 

a = (k 3 T n - 2£T 



k 3 T n - 2£T 

M-3/2 



From dM]) and ([1)5]), we get 

a k 

6 ~ 6 



k 3 nT 



-£l+n 



£(n - 2) 



(64) 
(65) 

(66) 



In subcase A = and £ = 0, metric (|57[) with expressions p, p, 9 and cr for this 
model are same as that of solution (f3~4"]l of Camci et al. [ID]. If we set £o = 0, 
this metric (|5"T|) represents the solution obtained by Pradhan et al. [H] . 

In subcase A = 0,£ = 0,n = 0, after an inverse transformation, metric (|57p 
reduces to the form 



(67) 



ds 2 = dt 2 - K 5 (t + U) 2/3 dx 2 - e 2ax \dy 2 + K%(t + U) 4/3 dz 2 



where t\ is an integrating constant. This model is different from the model (|36 
by a change of scale. 

In subcase A = 0, k = 0, same model as (|3T| is obtained. 



Further in subcase A = 0,£ = 0/c3 = 0, we see that the metric ([ST]) takes 
the form 

ds 2 = dt 2 - 2£K 6 (t + t b ) 2 ' 3 [dx 2 + e 2ax (bt m *dy 2 + b-H- m2 dz 2 )] , (68) 

where t 5 is a constant, bt m2 = 2£Kj^ ' (t + t 5 ) 2 &-$ and K 6 = 3{2 ^ . 

This metric is only different from p7|) by a change of sign. Also, in each of 
subcase the physical and kinematical properties of obtained metric are same as 
that of Case(I). Therefore, we do not consider here them. 
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3.4 Case (IV) : Let B = f n , where n is any real number. 

In this case Eq. (f2Tj) gives 



C = k 4 f n exp 
and then from (fT4"]l . we obtain 

A 2 = k 4 r 2n exp 
Hence the metric ((H) reduces to 



ds 2 - ^ 2n 



t- exp I — r 



f "exp( 



dy + exp 



2k 

V 



dz' 



dx z 



(69) 



(70) 



(71) 



where, without any loss of generality, the constant k 4 is taken equal to 1. Ex- 
pressions for physical and kinematical parameters for the model (|T1[) are not 
given here, but it is observed that the properties of the metric (J7TJ are same as 
that of the solution (|42|) , i.e. the Case (II). 



4 Concluding remarks 

In this paper we have described a new exact solutions of EFES for Bianchi type 
V spacetime with a bulk viscous fluid as the source of matter and cosmological 
term A varying with time. Using a generation technique followed by Camci 
et al. 40], it is shown that the Einstein's field equations are solvable for any 
arbitrary cosmic scale function. Starting from particular cosmic functions, new 
classes of spatially homogeneous and anisotropic cosmological models have been 
investigated for which the fluids are acceleration and rotation free but they do 
have expansion and shear. For a = in the metric (j4|), we obtained metrics 
as LRS Bianchi type I model (Hajj-Boutros [33], Hajj-Boutros and Sfeila [34] . 
Ram [35], Mazumder [35], Pradhan and Kumar [37]) and Pradhan et al. [H]. It 
is also seen that the solutions obtained by Camci et al. [40] , Ram [35] , Pradhan 
and Kumar [37] and Bianchi type V models studied by Pradhan et al. [41] are 
particular cases (except one) of our solutions. 

The cosmological constants in all models given in Section 3 are decreasing 
functions of time (except one, Model II of Case (I)) and they all approach a 
small positive value at late times which are supported by the results from re- 
cent supernova observations recently obtained by the High-z Supernova Team 
and Supernova Cosmological Project (Garnavich et al. [9], Perlmutter et al. 
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[ID] , Riess ef aZ. [II], Schmidt et al.\FZ\). 



The features of these new solutions are that from our wide range of choice 
of parameters for Model I of Cases (I), (II), and (III), it is apparent from fig- 
ures that energy density p{t) and dynamic cosmological term A(t) are decreas- 
ing function of time, remain positive and small at later stage during evolution. 
These two quantities remain finite and do not become zero at later stage. Hence 
this seem to be physically viable to explore physical mechanism for further detail 
study depending on relevance of the physical problem. These studies show that 
the bulk viscous effect is apparent on A(t). In many astrophysical situations 
the viscosity calculated using statistical consideration do not give satisfactory 
viscous number. Hence the specific thing about bulk viscosity will depend on 
the details of the viscous nature of the matter. We have explored Model II of 
Cases (I) , (II) and (III) . So we feel that it is not necessary to display behaviour of 
p(t) and A(t). We have studied these models and no anomaly has been observed. 

Model II of Case (I) is explored for variety of parameters. We find that 
many cases energy density decreases sharply and becomes negative very fast 
and later stage also remains negative but there is a slight increase at later stage 
and remains constant (negative constant), which we feel undesirable feature. So 
it is not discussed in the paper. Also in many cases it initially oscillatory. In this 
case A(i) shows very peculiar behaviour. Only in few cases A is initially negative 
then becomes positive but in all cases during initial evolution A has singularities 
near t close to zero and a few at later finite time. The cause of this behaviour 
is not very apparent from the model. Otherwise in wide range of parameters A 
shows eratic behaviour. Most of the cases singular and eratic behaviour of A(t) 
is persistent and cause is not known. So we do not discuss further. This may 
require detailed study separately and this paper is not appropriate place. 
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Figure 1: (Top) Plot of p — > t and (bottom) A — > r for parameters n 
fci = 1, 7 = 0.5, a = 1, £o = 1 and rest of the constants are set to 1. 
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Figure 2: (Top) Plot of p — ► f and (bottom) A — ► f for parameters n = 0.45, 
K = 2.0, fci = 2, 7 = 0.5, a = 1, £o = 1, M = 1 and rest of the constants are 
set to 1. 
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Figure 3: (Top) Plot of p — > T and (bottom) A — > T for parameters n = 0.45, 
K = 2.0, &3 = 1, 7 = 0.5, a = 1, £o = 1 an d rest of the constants are set to 1. 



22 



